Noise and diffusion of particles obeying asymmetric exclusion processes 
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The relation between noise and Pick's diffusion coefficient in barrier limited transport associated 
with hopping or tunneling mechanisms of particles obeying the asymmetric simple exclusion pro- 
cesses (ASEP) is physically assessed by Monte Carlo simulations. Por a closed ring consisting of a 
large number of barriers the diffusion coefficient is related explicitly to the current noise thus reveal- 
ing the existence of a generalized Nyquist-Einstein relation. Both diffusion and noise are confirmed 
to decrease as the square root of the number of barriers as a consequence of the correlation induced 
by ASEP. By contrast, for an open linear chain of barriers the diffusion coefficient is found to be no 
longer related to current noise. Here diffusion depends on particle concentration but is independent 
of the number of barriers. 
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The inter-relations between noise and diffusion in 
charge transport is a pillar of non-equilibrium kinetics 
since the theoretical interpretation of the Brownian mo- 
tion by Einsteini For a kinetic described within a con- 

\ tinuous transport model, where quasi-particles undergo 
local scattering events between stochastic free flights, 
the inter-relation between noise and diffusion was inves- 
tigated by a number of theoretical approaches ranging 
from analytical modela^i^i^i^ to numerical solutions of 

I the appropriate kinetic equationsi^ The case of a bar- 
rier transport model, dominated by tunnelling and/or 
hopping processes, is less developed. Here, noise was 
mostly investigated for the case of single and multiple 
quantum barriers i2ii24iiii2iiii4 By contrast, a few semi- 
nal works have tackled the problem of noise in hopping 
systemsiSii^ and that of diffusion in both tunneling and 
hopping systems liii^iiSiSS For the case of a very large 
number of barriers, the asymmetric simple exclusion pro- 
cess (ASEP) has been widely used in the recent past 
as a relevant physical model for the description of non- 
equilibrium dynamics i2ii22iSi2ii2^ In this context, two 
systems of basic interest are the closed ring and the open 
linear chain consisting of a set of multiple barriers, which 
are the prototypes of closed and open systems driven by 
hopping or tunneling transport mechanisms. Here, dif- 
fusion was investigated in by analytical meanSf2iiS& and 
current noise with Monte Carlo simulations. -'^^'^^ How- 
ever, the attempt to interrelate diffusion and current 
noise in these systems remains a largely unexplored is- 
sue. In particular, the dependence or less of diffusion on 
the number of barriers, as well as the prediction of a par- 
tial or complete suppression of shot noise in the presence 
of a large number of barriers are intriguing features still 
lacking of a microscopic interpretation. ^'"^'^^'^^ The aim 
of the present work is to address this issue by first prin- 
ciples Monte Carlo simulations. Accordingly, diffusion 
is obtained by the calculation of the time evolution of 
the spreading in space of a particle ensemble and current 
noise by the calculation of the autocorrelation function 
of current fluctuations as measured in the outside cir- 



cuit. The main features of diffusion and noise and their 
interrelation are thus quantitatively assessed on a kinetic 
physical ground. 

We take a physical system consisting of a number Nyj of 
hopping sites, separated by a constant distance I, whose 
total length is L = N^l. The tunneling rate between 
two adjacent site is assumed to have the same value F. 
We then consider alternatively periodic- (closed ring) or 
open-boundary (linear chain) conditions. By imposing a 
current determined by a transition rate F — 10^^ s~^ to 
a structure with Z = 3.2 nm (which are taken as plausible 
parameters for a real case) we evaluate current, diffusion, 
and noise making use of an ensemble Monte Carlo simu- 
lator. It is convenient to deflne the dimensionless carrier 
concentration p = {N)/Ni^, where (TV) is the average 
number of carrier inside the sample. We then introduce 
correlations between carriers within the ASEP model by 
imposing a maximum occupation number v for each site. 
In particular, when = 1 if a site is occupied by one 
carrier then no other carrier can jump to this site, thus 
carriers are strongly correlated. When v oo a site can 
be occupied by an arbitrary number of carriers, thus car- 
riers are totally uncorrelated. The instantaneous current 
I{t) is calculated as^ 

N(t) 

where e is the unit charge, N{t) the instantaneous num- 
ber of carriers inside the structure, Vi{t) the instanta- 
neous velocity of the i-th carrier, Vdit) the instantaneous 
drift velocity of the carrier ensemble. Under steady state 
I{t) is a stochastic variable that accounts for fluctuations 
in carrier number and velocity. In particular, for our dis- 
crete system Vi{t) = lS£,i/St, with the position index 
of the i-th particleiii 

For the (longitudinal) diffusion coefficient D, follow- 
ing Fick's law we make use of its definition as a spatial 
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FIG. 1: (Color online) Cloded ring. Comparison between 
the analytical diffusion coefficient in Eqs. Q and Q (dashed 
lines) and that obtained from simulations (symbols) 
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FIG. 2: (Color online) Closed ring. Comparison between the 
analytical Fano factor in Eqs. Q, IIH (dashed lines) and 
that obtained from simulations (symbols). 
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where Az{t) is the distance between the final and the 
initial hopping sites, brackets mean average over a sta- 
tistical ensemble (up to 10^) of identical systems, and the 
time derivative is carried out in a time domain shorter 
than the transit time but sufficiently long to extrapolate 
the long time limit for which D is found to be indepen- 
dent of time. 

The spectral density of current fiuctuations at zero fre- 
quency is2& 



5/ = 4 / dt{siio)si{t)) = s";" + + s';"" + s^j 

Jo 

(3) 

where Sf' and {S^'"" + S^/^) refer to the three con- 
tributions (drift velocity, number and cross-correlations 
between them) in which the total spectral density can 
be decomposed. With Monte Carlo simulations these 
terms can be calculated separately. The Fano factor is 
7 = S,/(2e(7)). 

A. Closed ring. The periodic boundary conditions ad- 
equate to this structure consists in imposing that, for a 



finite number of sites, the last site is directly connected to 
the first one by the same hopping rate F. Let us consider 
the case of correlated carriers (i^ = 1). For a given car- 
rier concentration and for large N.^ analytical theorjii^iSi 
gives the following predictions: for the average current 



(I),=eTpil-p) , 



(4) 



and for the diffusion coefficient Ai (following ^Tj, with 
the subscript 1 labeling the case of the ring geometry). 



Ai 
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(5) 



[We suppose that in general the value of Ai differs from 
that of D in Eq.©.] 

The results of the simulations for the diffusion coeffi- 
cient are reported in Fig. ^ Here, the identity Ai — D 
is confirmed for all the concentrations considered. 

For the current noise (in this case due only to veloc- 
ity fluctuations since the number of carriers is rigorously 
constant in time) the simulations confirm the relation 



Si = S'"/ = ^p'Ai 



4e^ 



(6) 



with the corresponding Fano factor (Fig. O, given by 
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Equation © , by revealing a strict relation between noise 
and diffusion, takes the form of a generalized Nyquist- 
Einstein relation2i2£. The reason of diffusion and noise 
suppression as l/y/N^, confirmed by the simulations, is 
attributed to the strong correlation among carriers. To 
support this interpretation, we considered also the case of 
uncorrelated carriers (i.e. in the absence of ASEP) where, 
for a given carrier concentration, analytical theory gives 
the following predictions: for the average currenlji^ 



(/)o = eVp , 
for the diffusion coefficientiSiiS 

Dn = — . 



(8) 



(9) 



(with the subscript labeling the case of uncorrelated 
particles) and for the current noise the standard Nyquist 
Einstein relation^S 



Si = S^/ = 



4e2 ^Do 



with the corresponding Fano factor 
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1 



(10) 



(11) 



The result of simulations confirms that, in the absence 
of ASEP, diffusion becomes independent of iV^, (see the 
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FIG. 3: (Color online) Linear chain. Spreading diffusion co- 
efficient obtained from simulations in the presence of ASEP 
and for different input and output rates. The value for uncor- 
related carriers is reported for comparison. Curves are guide 
to the eyes. 
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FIG. 4: (Color online) Linear chain. Comparison between 
analytical Fano factors in Eq. 11161 and full shot noise (dashed 
lines), with that obtained from simulations (symbols). 



curve uncorrelated in Fig. and the Fano factor de- 
creases as l/Nw (see the curve uncorrelated in Fig. [2J). 
In all the cases considered here (even "when ~ 1 so 
that the non-passing constraint is accompHshcd) the time 
evolution of the variance in space of the carrier ensem- 
ble is found to be linear, contrarily to the suggestion of 
a sub-diffusive (and thus sublinear) behaviori2i On the 
other hand, the suppression of diffusion in the presence of 
ASEP is found to be related to a subpoissian variance of 
the number of scattering events, "which is found to exhibit 
the expected behavior. 

B. Open linear chain. The boundary conditions ad- 
equate to this structure consists in connecting the t"wo 
terminals of the device "with t"wo reservoirs, "where Tin x T 
and Tout x T are the transition rates from the left reser- 
voir to the first site of the device, and from the last site 
of the device to the right reservoir, respectively. For con- 
venience, the values of Tin(^out) are taken in the range 
bet"ween and 1 bein g equi valent respectively to the a 
and f3 parameters in |2ll2(il | and to /l and (1 — in 

m 

Let us first consider the case of correlated carriers 
(ASEP model). For a given carrier concentration and 



for large N^, analytical theory gives for the average cur- 
rent the same expression of Eq. Q). The diffusion co- 
efficient of the linear chain A2 (follo"wing [2G|, "with the 
subscript 2 labeling this specific case) takes the form: if 
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[Again A2 is supposed to differ in general from D ob- 
tained from Eq.®.] The different analytical expressions 
for A2 are a consequence of the different values taken by 
^in,out and, in turn, by p in the steady state. Indeed, 
the tuning of the Ti^^out controls the strength of the cor- 
relation among carriers induced by ASEP and thus the 
particle density p inside the device, as summarized in the 
phase diagram reported in Fig. 2 of [2l| . 

Figure reports the Pick's diffusion coefficient D ob- 
tained from simulations for the case of the linear chain. 
Here, D is found to be practically independent of Nuj- 
Furthermore, in the presence of ASEP the value of D is 
systematically lo"wer than that of uncorrelated particles 
Do ■ For the case Tin + Tout = 1 , the value of the diffusion 
coefficient is well described by the relation 



D 



PP, 



PT 



(1-p) 



(14) 



We notice that in the above expression, D becomes van- 
ishing small for Tout ~* because in this limit spreading 
and current of carriers through the structure tends to 
stop. In the presence of ASEP the values of diffusion 
obtained by simulations are found to differ significantly 
from those given by analytical expressions, thus imply- 
ing that the quantities A2 an D describe different micro- 
scopic processes. 

Figure ^ reports the results of the simulations for the 
current noise [in this case due to the sum of all the con- 
tributions in Eq. ©J in terms of the Fano factor. From 
simulations, "within numerical uncertainty we find: 



Si = S"/ + Sf + S'/"^ + S^j 



4e2 ^ 



(15) 



with the corresponding Fano factor satisfying the relation 

2A2 
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From the above expressions we conclude that A2 de- 
scribes the total current noise instead of the Pick's dif- 
fusion process. The quantity A2 is found to agree well 
with available analytical expressions, as predicted by Eq. 
(|15|l , in full agreement with the results of (2^ . The quan- 
tity D is found to depend upon the degree of correla- 
tion, to be independent of the number of sites and, when 
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FIG. 5: (Color online) Different contributions and total noise 
power for the case of ASEP in the open linear with Fin = 
Tout ~ 1- Curves are guide to the eyes 

Tin + ToMt = 1, to be proportional to Tout- Since A2 is 
related to the number of particles that entered the device 
up to time we conclude that this stochastic quantity 
accounts for both velocity and number fluctuations. 

By turning off the ASEP correlation, for a given carrier 
concentration, analytical theories give: for the average 
current Eq. (jS]), for the diffusion coefficient Eq. and 
for the current noise (full shot noise is expected): 

Sj^Sf ^^{vd)'{N)TTR^2e{I) (17) 

which follows from a correlation function with triangular 
shape vanishing at the transit time ttr — L/{vd). As 
a consequence 7 = 1, as confirmed by the results of the 
simulations reported in Fig. ^ 

Finally, we have investigated separately the contribu- 
tions to the total current noise, which come from velocity, 
number, and their cross-correlations in the presence of 
ASEP. From the results © we can see that different con- 



tributions are comparable in magnitude, with the cross 
terms, which are responsible of shot noise suppression, 
being negative. We further notice that the comparison 
between the two uncorrelated current-noise levels belong- 
ing to the closed ring (Fig. 12) and to the open linear 
chain (Fig. ^ shows that shot noise of the linear chain 
exceeds velocity noise of the closed ring for the ratio L/l, 
as expected. 

In conclusion, we have carried out a simulative investi- 
gation of the inter-relations between noise and diffusion 
in barrier limited transport under the ASEP condition. 
For the case of the closed ring, since the number of parti- 
cles is fixed, only the noise related to velocity fluctuations 
is present. Here, the diffusion coefficient obtained from 
the Pick's law is explicitly related to current noise, both 
in the presence and in the absence of the ASEP. There- 
fore, evidence is provided for the existence of a gener- 
alized Nyquist-Einstein relation allowing the determina- 
tion of diffusion from a noise measurement or viceversa. 
The correlations introduced by ASEP are found to be 
responsible of the dependence of diffusion upon the in- 
verse square root of the device length. For the case of 
the open linear chain the diffusion coefficient obtained 
from Pick's law is not related to the current noise, which 
now contains contributions coming from velocity, num- 
ber and their cross-correlations. Here the diffusion coef- 
ficient is found to be independent of the number of sites 
but to depend on the strength of the correlation that is 
ultimately controlled by the carrier density. We remark, 
that the total current noise is found to be related to a 
"diffusion constant" (more properly a counting statistics 
property^), whose analytical expressions are satisfacto- 
rily confirmed by simulations. 
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